MHC class II genes mediate susceptibility and resistance to coronavirus infections in bats.

Understanding the immunogenetic basis of coronavirus (CoV) susceptibility in major pathogen reservoirs, such as bats, is central to inferring their zoonotic potential. Members of the cryptic Hipposideros bat species complex differ in CoV susceptibility, but the underlying mechanisms remain unclear. The genes of the major histocompatibility complex (MHC) are the best understood genetic basis of pathogen resistance, and differences in MHC diversity are one possible reason for asymmetrical infection patterns among closely related species. Here, we aimed to link asymmetries in observed CoV (CoV-229E, CoV-2B and CoV-2Bbasal) susceptibility to immunogenetic differences amongst four Hipposideros bat species. From the 2072 bats assigned to their respective species using the mtDNA cytochrome b gene, members of the most numerous and ubiquitous species, Hipposideros caffer D, were most infected with CoV-229E and SARS-related CoV-2B. Using a subset of 569 bats, we determined that much of the existent allelic and functional (i.e. supertype) MHC DRB class II diversity originated from common ancestry. One MHC supertype shared amongst all species, ST12, was consistently linked to susceptibility with CoV-229E, which is closely related to the common cold agent HCoV-229E, and infected bats and those carrying ST12 had a lower body condition. The same MHC supertype was connected to resistance to CoV-2B, and bats with ST12 were less likely be co-infected with CoV-229E and CoV-2B. Our work suggests a role of immunogenetics in determining CoV susceptibility in bats. We advocate for the preservation of functional genetic and species diversity in reservoirs as a means of mitigating the risk of disease spillover.

Next Generation Sequencing Revolutionizes Organismal Biology Research in Bats.

The advent of next generation sequencing technologies (NGS) has greatly accelerated our understanding of critical aspects of organismal biology from non-model organisms. Bats form a particularly interesting group in this regard, as genomic data have helped unearth a vast spectrum of idiosyncrasies in bat genomes associated with bat biology, physiology, and evolution. Bats are important bioindicators and are keystone species to many eco-systems. They often live in proximity to humans and are frequently associated with emerging infectious diseases, including the COVID-19 pandemic. Nearly four dozen bat genomes have been published to date, ranging from drafts to chromosomal level assemblies. Genomic investigations in bats have also become critical towards our understanding of disease biology and host-pathogen coevolution. In addition to whole genome sequencing, low coverage genomic data like reduced representation libraries, resequencing data, etc. have contributed significantly towards our understanding of the evolution of natural populations, and their responses to climatic and anthropogenic perturbations. In this review, we discuss how genomic data have enhanced our understanding of physiological adaptations in bats (particularly related to ageing, immunity, diet, etc.), pathogen discovery, and host pathogen co-evolution. In comparison, the application of NGS towards population genomics, conservation, biodiversity assessment, and functional genomics has been appreciably slower. We reviewed the current areas of focus, identifying emerging topical research directions and providing a roadmap for future genomic studies in bats.

Bat Red Blood Cells Express Nucleic Acid-Sensing Receptors and Bind RNA and DNA.

RBCs demonstrate immunomodulatory capabilities through the expression of nucleic acid sensors. However, little is known about bat RBCs, and no studies have examined the immune function of bat erythrocytes. In this study, we show that bat RBCs express the nucleic acid-sensing TLRs TLR7 and TLR9 and bind the nucleic acid ligands, ssRNA, and CpG DNA. Collectively, these data suggest that, like human RBCs, bat erythrocytes possess immune function and may be reservoirs for nucleic acids. These findings provide unique insight into bat immunity and may uncover potential mechanisms by which virulent pathogens of humans are concealed in bats.

Analogous comparison unravels heightened antiviral defense and boosted viral infection upon immunosuppression in bat organoids.

Horseshoe bats host numerous SARS-related coronaviruses without overt disease signs. Bat intestinal organoids, a unique model of bat intestinal epithelium, allow direct comparison with human intestinal organoids. We sought to unravel the cellular mechanism(s) underlying bat tolerance of coronaviruses by comparing the innate immunity in bat and human organoids. We optimized the culture medium, which enabled a consecutive passage of bat intestinal organoids for over one year. Basal expression levels of IFNs and IFN-stimulated genes were higher in bat organoids than in their human counterparts. Notably, bat organoids mounted a more rapid, robust and prolonged antiviral defense than human organoids upon Poly(I:C) stimulation. TLR3 and RLR might be the conserved pathways mediating antiviral response in bat and human intestinal organoids. The susceptibility of bat organoids to a bat coronavirus CoV-HKU4, but resistance to EV-71, an enterovirus of exclusive human origin, indicated that bat organoids adequately recapitulated the authentic susceptibility of bats to certain viruses. Importantly, TLR3/RLR inhibition in bat organoids significantly boosted viral growth in the early phase after SARS-CoV-2 or CoV-HKU4 infection. Collectively, the higher basal expression of antiviral genes, especially more rapid and robust induction of innate immune response, empowered bat cells to curtail virus propagation in the early phase of infection.

Insertion/deletion hotspots in the Nsp2, Nsp3, S1, and ORF8 genes of SARS-related coronaviruses.

The genome of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) contains many insertions/deletions (indels) from the genomes of other SARS-related coronaviruses. Some of the identified indels have recently reported to involve relatively long segments of 10-300 consecutive bases and with diverse RNA sequences around gaps between virus species, both of which are different characteristics from the classical shorter in-frame indels. These non-classical complex indels have been identified in non-structural protein 3 (Nsp3), the S1 domain of the spike (S), and open reading frame 8 (ORF8). To determine whether the occurrence of these non-classical indels in specific genomic regions is ubiquitous among broad species of SARS-related coronaviruses in different animal hosts, the present study compared SARS-related coronaviruses from humans (SARS-CoV and SARS-CoV-2), bats (RaTG13 and Rc-o319), and pangolins (GX-P4L), by performing multiple sequence alignment. As a result, indel hotspots with diverse RNA sequences of different lengths between the viruses were confirmed in the Nsp2 gene (approximately 2500-2600 base positions in the overall 29,900 bases), Nsp3 gene (approximately 3000-3300 and 3800-3900 base positions), N-terminal domain of the spike protein (21,500-22,500 base positions), and ORF8 gene (27,800-28,200 base positions). Abnormally high rate of point mutations and complex indels in these regions suggest that the occurrence of mutations in these hotspots may be selectively neutral or even benefit the survival of the viruses. The presence of such indel hotspots has not been reported in different human SARS-CoV-2 strains in the last 2 years, suggesting a lower rate of indels in human SARS-CoV-2. Future studies to elucidate the mechanisms enabling the frequent development of long and complex indels in specific genomic regions of SARS-related coronaviruses would offer deeper insights into the process of viral evolution.

RNA structure-altering mutations underlying positive selection on Spike protein reveal novel putative signatures to trace crossing host-species barriers in Betacoronavirus.

Similar to other RNA viruses, the emergence of Betacoronavirus relies on cross-species viral transmission, which requires careful health surveillance monitoring of protein-coding information as well as genome-wide analysis. Although the evolutionary jump from natural reservoirs to humans may be mainly traced-back by studying the effect that hotspot mutations have on viral proteins, it is largely unexplored if other impacts might emerge on the structured RNA genome of Betacoronavirus. In this survey, the protein-coding and viral genome architecture were simultaneously studied to uncover novel insights into cross-species horizontal transmission events. We analysed 1,252,952 viral genomes of SARS-CoV, MERS-CoV, and SARS-CoV-2 distributed across the world in bats, intermediate animals, and humans to build a new landscape of changes in the RNA viral genome. Phylogenetic analyses suggest that bat viruses are the most closely related to the time of most recent common ancestor of Betacoronavirus, and missense mutations in viral proteins, mainly in the S protein S1 subunit: SARS-CoV (G > T; A577S); MERS-CoV (C > T; S746R and C > T; N762A); and SARS-CoV-2 (A > G; D614G) appear to have driven viral diversification. We also found that codon sites under positive selection on S protein overlap with non-compensatory mutations that disrupt secondary RNA structures in the RNA genome complement. These findings provide pivotal factors that might be underlying the eventual jumping the species barrier from bats to intermediate hosts. Lastly, we discovered that nearly half of the Betacoronavirus genomes carry highly conserved RNA structures, and more than 90% of these RNA structures show negative selection signals, suggesting essential functions in the biology of Betacoronavirus that have not been investigated to date. Further research is needed on negatively selected RNA structures to scan for emerging functions like the potential of coding virus-derived small RNAs and to develop new candidate antiviral therapeutic strategies.

Exceptional diversity and selection pressure on coronavirus host receptors in bats compared to other mammals.

Pandemics originating from non-human animals highlight the need to understand how natural hosts have evolved in response to emerging human pathogens and which groups may be susceptible to infection and/or potential reservoirs to mitigate public health and conservation concerns. Multiple zoonotic coronaviruses, such as severe acute respiratory syndrome-associated coronavirus (SARS-CoV), SARS-CoV-2 and Middle Eastern respiratory syndrome-associated coronavirus (MERS-CoV), are hypothesized to have evolved in bats. We investigate angiotensin-converting enzyme 2 (ACE2), the host protein bound by SARS-CoV and SARS-CoV-2, and dipeptidyl-peptidase 4 (DPP4 or CD26), the host protein bound by MERS-CoV, in the largest bat datasets to date. Both the ACE2 and DPP4 genes are under strong selection pressure in bats, more so than in other mammals, and in residues that contact viruses. Additionally, mammalian groups vary in their similarity to humans in residues that contact SARS-CoV, SARS-CoV-2 and MERS-CoV, and increased similarity to humans in binding residues is broadly predictive of susceptibility to SARS-CoV-2. This work augments our understanding of the relationship between coronaviruses and mammals, particularly bats, provides taxonomically diverse data for studies of how host proteins are bound by coronaviruses and can inform surveillance, conservation and public health efforts.

Genomic determinants of Furin cleavage in diverse European SARS-related bat coronaviruses.

The furin cleavage site (FCS) in SARS-CoV-2 is unique within the Severe acute respiratory syndrome-related coronavirus (SrC) species. We re-assessed diverse SrC from European horseshoe bats and analyzed the spike-encoding genomic region harboring the FCS in SARS-CoV-2. We reveal molecular features in SrC such as purine richness and RNA secondary structures that resemble those required for FCS acquisition in avian influenza viruses. We discuss the potential acquisition of FCS through molecular mechanisms such as nucleotide substitution, insertion, or recombination, and show that a single nucleotide exchange in two European bat-associated SrC may suffice to enable furin cleavage. Furthermore, we show that FCS occurrence is variable in bat- and rodent-borne counterparts of human coronaviruses. Our results suggest that furin cleavage sites can be acquired in SrC via conserved molecular mechanisms known in other reservoir-bound RNA viruses and thus support a natural origin of SARS-CoV-2.

Variation in synonymous nucleotide composition among genomes of sarbecoviruses and consequences for the origin of COVID-19.

The subgenus Sarbecovirus includes two human viruses, SARS-CoV and SARS-CoV-2, respectively responsible for the SARS epidemic and COVID-19 pandemic, as well as many bat viruses and two pangolin viruses. Here, the synonymous nucleotide composition (SNC) of Sarbecovirus genomes was analysed by examining third codon-positions, dinucleotides, and degenerate codons. The results show evidence for the eight following groups: (i) SARS-CoV related coronaviruses (SCoVrC including many bat viruses from China), (ii) SARS-CoV-2 related coronaviruses (SCoV2rC; including five bat viruses from Cambodia, Thailand and Yunnan), (iii) pangolin sarbecoviruses, (iv) three bat sarbecoviruses showing evidence of recombination between SCoVrC and SCoV2rC genomes, (v) two highly divergent bat sarbecoviruses from Yunnan, (vi) the bat sarbecovirus from Japan, (vii) the bat sarbecovirus from Bulgaria, and (viii) the bat sarbecovirus from Kenya. All these groups can be diagnosed by specific nucleotide compositional features except the one concerned by recombination between SCoVrC and SCoV2rC. In particular, SCoV2rC genomes have less cytosines and more uracils at third codon-positions than other sarbecoviruses, whereas the genomes of pangolin sarbecoviruses show more adenines at third codon-positions. I suggest that taxonomic differences in the imbalanced nucleotide pools available in host cells during viral replication can explain the eight groups of SNC here detected among Sarbecovirus genomes. A related effect due to hibernating bats and their latitudinal distribution is also discussed. I conclude that the two independent host switches from Rhinolophus bats to pangolins resulted in convergent mutational constraints and that SARS-CoV-2 emerged directly from a horseshoe bat sarbecovirus.

Insertion-and-Deletion Mutations between the Genomes of SARS-CoV, SARS-CoV-2, and Bat Coronavirus RaTG13.

The evolutional process of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) development remains inconclusive. This study compared the genome sequences of severe acute respiratory syndrome coronavirus (SARS-CoV), bat coronavirus RaTG13, and SARS-CoV-2. In total, the genomes of SARS-CoV-2 and RaTG13 were 77.9% and 77.7% identical to the genome of SARS-CoV, respectively. A total of 3.6% (1,068 bases) of the SARS-CoV-2 genome was derived from insertion and/or deletion (indel) mutations, and 18.6% (5,548 bases) was from point mutations from the genome of SARS-CoV. At least 35 indel sites were confirmed in the genome of SARS-CoV-2, in which 17 were with ≥10 consecutive bases long. Ten of these relatively long indels were located in the spike (S) gene, five in nonstructural protein 3 (Nsp3) gene of open reading frame (ORF) 1a, and one in ORF8 and noncoding region. Seventeen (48.6%) of the 35 indels were based on insertion-and-deletion mutations with exchanged gene sequences of 7-325 consecutive bases. Almost the complete ORF8 gene was replaced by a single 325 consecutive base-long indel. The distribution of these indels was roughly in accordance with the distribution of the rate of point mutation rate around the indels. The genome sequence of SARS-CoV-2 was 96.0% identical to that of RaTG13. There was no long insertion-and-deletion mutation between the genomes of RaTG13 and SARS-CoV-2. The findings of the uneven distribution of multiple indels and the presence of multiple long insertion-and-deletion mutations with exchanged consecutive base sequences in the viral genome may provide insights into SARS-CoV-2 development. IMPORTANCE The developmental mechanism of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) remains inconclusive. This study compared the base sequence one-by-one between severe acute respiratory syndrome coronavirus (SARS-CoV) or bat coronavirus RaTG13 and SARS-CoV-2. The genomes of SARS-CoV-2 and RaTG13 were 77.9% and 77.7% identical to the genome of SARS-CoV, respectively. Seventeen of the 35 sites with insertion and/or deletion mutations between SARS-CoV-2 and SARS-CoV were based on insertion-and-deletion mutations with the replacement of 7-325 consecutive bases. Most of these long insertion-and-deletion sites were concentrated in the nonstructural protein 3 (Nsp3) gene of open reading frame (ORF) 1a, S1 domain of the spike protein, and ORF8 genes. Such long insertion-and-deletion mutations were not observed between the genomes of RaTG13 and SARS-CoV-2. The presence of multiple long insertion-and-deletion mutations in the genome of SARS-CoV-2 and their uneven distributions may provide further insights into the development of the virus.

Crystal Structures of Bat and Human Coronavirus ORF8 Protein Ig-Like Domain Provide Insights Into the Diversity of Immune Responses.

ORF8 is a viral immunoglobulin-like (Ig-like) domain protein encoded by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) RNA genome. It tends to evolve rapidly and interfere with immune responses. However, the structural characteristics of various coronavirus ORF8 proteins and their subsequent effects on biological functions remain unclear. Herein, we determined the crystal structures of SARS-CoV-2 ORF8 (S84) (one of the epidemic isoforms) and the bat coronavirus RaTG13 ORF8 variant at 1.62 Å and 1.76 Å resolution, respectively. Comparison of these ORF8 proteins demonstrates that the 62-77 residues in Ig-like domain of coronavirus ORF8 adopt different conformations. Combined with mutagenesis assays, the residue Cys20 of ORF8 is responsible for forming the covalent disulfide-linked dimer in crystal packing and in vitro biochemical conditions. Furthermore, immune cell-binding assays indicate that various ORF8 (SARS-CoV-2 ORF8 (L84), ORF8 (S84), and RaTG13 ORF8) proteins have different interaction capabilities with human CD14+ monocytes in human peripheral blood. These results provide new insights into the specific characteristics of various coronavirus ORF8 and suggest that ORF8 variants may influence disease-related immune responses.

Genetic identification of bat species for pathogen surveillance across France.

With more than 1400 chiropteran species identified to date, bats comprise one-fifth of all mammalian species worldwide. Many studies have associated viral zoonoses with 45 different species of bats in the EU, which cluster within 5 families of bats. For example, the Serotine bats are infected by European Bat 1 Lyssavirus throughout Europe while Myotis bats are shown infected by coronavirus, herpesvirus and paramyxovirus. Correct host species identification is important to increase our knowledge of the ecology and evolutionary pattern of bat viruses in the EU. Bat species identification is commonly determined using morphological keys. Morphological determination of bat species from bat carcasses can be limited in some cases, due to the state of decomposition or nearly indistinguishable morphological features in juvenile bats and can lead to misidentifications. The overall objective of our study was to identify insectivorous bat species using molecular biology tools with the amplification of the partial cytochrome b gene of mitochondrial DNA. Two types of samples were tested in this study, bat wing punches and bat faeces. A total of 163 bat wing punches representing 22 species, and 31 faecal pellets representing 7 species were included in the study. From the 163 bat wing punches tested, a total of 159 were genetically identified from amplification of the partial cyt b gene. All 31 faecal pellets were genetically identified based on the cyt b gene. A comparison between morphological and genetic determination showed 21 misidentifications from the 163 wing punches, representing ~12.5% of misidentifications of morphological determination compared with the genetic method, across 11 species. In addition, genetic determination allowed the identification of 24 out of 25 morphologically non-determined bat samples. Our findings demonstrate the importance of a genetic approach as an efficient and reliable method to identify bat species precisely.

Constitutive IFNα Protein Production in Bats.

Bats are the only mammals with self-powered flight and account for 20% of all extant mammalian diversity. In addition, they harbor many emerging and reemerging viruses, including multiple coronaviruses, several of which are highly pathogenic in other mammals, but cause no disease in bats. How this symbiotic relationship between bats and viruses exists is not yet fully understood. Existing evidence supports a specific role for the innate immune system, in particular type I interferon (IFN) responses, a major component of antiviral immunity. Previous studies in bats have shown that components of the IFN pathway are constitutively activated at the transcriptional level. In this study, we tested the hypothesis that the type I IFN response in bats is also constitutively activated at the protein level. For this, we utilized highly sensitive Single Molecule (Simoa) digital ELISA assays, previously developed for humans that we adapted to bat samples. We prospectively sampled four non-native chiroptera species from French zoos. We identified a constitutive expression of IFNα protein in the circulation of healthy bats, and concentrations that are physiologically active in humans. Expression levels differed according to the species examined, but were not associated with age, sex, or health status suggesting constitutive IFNα protein expression independent of disease. These results confirm a unique IFN response in bat species that may explain their ability to coexist with multiple viruses in the absence of pathology. These results may help to manage potential zoonotic viral reservoirs and potentially identify new anti-viral strategies.

Utilizing the VirIdAl Pipeline to Search for Viruses in the Metagenomic Data of Bat Samples.

According to various estimates, only a small percentage of existing viruses have been discovered, naturally much less being represented in the genomic databases. High-throughput sequencing technologies develop rapidly, empowering large-scale screening of various biological samples for the presence of pathogen-associated nucleotide sequences, but many organisms are yet to be attributed specific loci for identification. This problem particularly impedes viral screening, due to vast heterogeneity in viral genomes. In this paper, we present a new bioinformatic pipeline, VirIdAl, for detecting and identifying viral pathogens in sequencing data. We also demonstrate the utility of the new software by applying it to viral screening of the feces of bats collected in the Moscow region, which revealed a significant variety of viruses associated with bats, insects, plants, and protozoa. The presence of alpha and beta coronavirus reads, including the MERS-like bat virus, deserves a special mention, as it once again indicates that bats are indeed reservoirs for many viral pathogens. In addition, it was shown that alignment-based methods were unable to identify the taxon for a large proportion of reads, and we additionally applied other approaches, showing that they can further reveal the presence of viral agents in sequencing data. However, the incompleteness of viral databases remains a significant problem in the studies of viral diversity, and therefore necessitates the use of combined approaches, including those based on machine learning methods.

Genomic Characterization of Diverse Bat Coronavirus HKU10 in Hipposideros Bats.

Bats have been identified as natural reservoirs of a variety of coronaviruses. They harbor at least 19 of the 33 defined species of alpha- and betacoronaviruses. Previously, the bat coronavirus HKU10 was found in two bat species of different suborders, Rousettus leschenaultia and Hipposideros pomona, in south China. However, its geographic distribution and evolution history are not fully investigated. Here, we screened this viral species by a nested reverse transcriptase PCR in our archived samples collected over 10 years from 25 provinces of China and one province of Laos. From 8004 bat fecal samples, 26 were found to be positive for bat coronavirus HKU10 (BtCoV HKU10). New habitats of BtCoV HKU10 were found in the Yunnan, Guangxi, and Hainan Provinces of China, and Louang Namtha Province in Laos. In addition to H. pomona, BtCoV HKU10 variants were found circulating in Aselliscus stoliczkanus and Hipposideros larvatus. We sequenced full-length genomes of 17 newly discovered BtCoV HKU10 strains and compared them with previously published sequences. Our results revealed a much higher genetic diversity of BtCoV HKU10, particularly in spike genes and accessory genes. Besides the two previously reported lineages, we found six novel lineages in their new habitats, three of which were located in Yunnan province. The genotypes of these viruses are closely related to sampling locations based on polyproteins, and correlated to bat species based on spike genes. Combining phylogenetic analysis, selective pressure, and molecular-clock calculation, we demonstrated that Yunnan bats harbor a gene pool of BtCoV HKU10, with H. pomona as a natural reservoir. The cell tropism test using spike-pseudotyped lentivirus system showed that BtCoV HKU10 could enter cells from human and bat, suggesting a potential interspecies spillover. Continuous studies on these bat coronaviruses will expand our understanding of the evolution and genetic diversity of coronaviruses, and provide a prewarning of potential zoonotic diseases from bats.

A prenylated dsRNA sensor protects against severe COVID-19.

Inherited genetic factors can influence the severity of COVID-19, but the molecular explanation underpinning a genetic association is often unclear. Intracellular antiviral defenses can inhibit the replication of viruses and reduce disease severity. To better understand the antiviral defenses relevant to COVID-19, we used interferon-stimulated gene (ISG) expression screening to reveal that 2'-5'-oligoadenylate synthetase 1 (OAS1), through ribonuclease L, potently inhibits severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). We show that a common splice-acceptor single-nucleotide polymorphism (Rs10774671) governs whether patients express prenylated OAS1 isoforms that are membrane-associated and sense-specific regions of SARS-CoV-2 RNAs or if they only express cytosolic, nonprenylated OAS1 that does not efficiently detect SARS-CoV-2. In hospitalized patients, expression of prenylated OAS1 was associated with protection from severe COVID-19, suggesting that this antiviral defense is a major component of a protective antiviral response.

Orthogonal genome-wide screens of bat cells identify MTHFD1 as a target of broad antiviral therapy.

Bats are responsible for the zoonotic transmission of several major viral diseases, including those leading to the 2003 SARS outbreak and likely the ongoing COVID-19 pandemic. While comparative genomics studies have revealed characteristic adaptations of the bat innate immune system, functional genomic studies are urgently needed to provide a foundation for the molecular dissection of the viral tolerance in bats. Here we report the establishment of genome-wide RNA interference (RNAi) and CRISPR libraries for the screening of the model megabat, Pteropus alecto. We used the complementary RNAi and CRISPR libraries to interrogate P. alecto cells for infection with two different viruses: mumps virus and influenza A virus, respectively. Independent screening results converged on the endocytosis pathway and the protein secretory pathway as required for both viral infections. Additionally, we revealed a general dependence of the C1-tetrahydrofolate synthase gene, MTHFD1, for viral replication in bat cells and human cells. The MTHFD1 inhibitor, carolacton, potently blocked replication of several RNA viruses, including SARS-CoV-2. We also discovered that bats have lower expression levels of MTHFD1 than humans. Our studies provide a resource for systematic inquiry into the genetic underpinnings of bat biology and a potential target for developing broad-spectrum antiviral therapy.

Coronavirus surveillance in wildlife from two Congo basin countries detects RNA of multiple species circulating in bats and rodents.

Coronaviruses play an important role as pathogens of humans and animals, and the emergence of epidemics like SARS, MERS and COVID-19 is closely linked to zoonotic transmission events primarily from wild animals. Bats have been found to be an important source of coronaviruses with some of them having the potential to infect humans, with other animals serving as intermediate or alternate hosts or reservoirs. Host diversity may be an important contributor to viral diversity and thus the potential for zoonotic events. To date, limited research has been done in Africa on this topic, in particular in the Congo Basin despite frequent contact between humans and wildlife in this region. We sampled and, using consensus coronavirus PCR-primers, tested 3,561 wild animals for coronavirus RNA. The focus was on bats (38%), rodents (38%), and primates (23%) that posed an elevated risk for contact with people, and we found coronavirus RNA in 121 animals, of which all but two were bats. Depending on the taxonomic family, bats were significantly more likely to be coronavirus RNA-positive when sampled either in the wet (Pteropodidae and Rhinolophidae) or dry season (Hipposideridae, Miniopteridae, Molossidae, and Vespertilionidae). The detected RNA sequences correspond to 15 alpha- and 6 betacoronaviruses, with some of them being very similar (>95% nucleotide identities) to known coronaviruses and others being more unique and potentially representing novel viruses. In seven of the bats, we detected RNA most closely related to sequences of the human common cold coronaviruses 229E or NL63 (>80% nucleotide identities). The findings highlight the potential for coronavirus spillover, especially in regions with a high diversity of bats and close human contact, and reinforces the need for ongoing surveillance.

The Evolution of Severe Acute Respiratory Syndrome Coronavirus-2 during Pandemic and Adaptation to the Host.

Severe Acute Respiratory Syndrome Coronavirus-2 is a zoonotic virus with a possible origin in bats and potential transmission to humans through an intermediate host. When zoonotic viruses jump to a new host, they undergo both mutational and natural selective pressures that result in non-synonymous and synonymous adaptive changes, necessary for efficient replication and rapid spread of diseases in new host species. The nucleotide composition and codon usage pattern of SARS-CoV-2 indicate the presence of a highly conserved, gene-specific codon usage bias. The codon usage pattern of SARS-CoV-2 is mostly antagonistic to human and bat codon usage. SARS-CoV-2 codon usage bias is mainly shaped by the natural selection, while mutational pressure plays a minor role. The time-series analysis of SARS-CoV-2 genome indicates that the virus is slowly evolving. Virus isolates from later stages of the outbreak have more biased codon usage and nucleotide composition than virus isolates from early stages of the outbreak.

Mutations derived from horseshoe bat ACE2 orthologs enhance ACE2-Fc neutralization of SARS-CoV-2.

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike (S) protein mediates infection of cells expressing angiotensin-converting enzyme 2 (ACE2). ACE2 is also the viral receptor of SARS-CoV (SARS-CoV-1), a related coronavirus that emerged in 2002-2003. Horseshoe bats (genus Rhinolophus) are presumed to be the original reservoir of both viruses, and a SARS-like coronavirus, RaTG13, closely related to SARS-CoV-2, has been identified in one horseshoe-bat species. Here we characterize the ability of the S-protein receptor-binding domains (RBDs) of SARS-CoV-1, SARS-CoV-2, pangolin coronavirus (PgCoV), RaTG13, and LyRa11, a bat virus similar to SARS-CoV-1, to bind a range of ACE2 orthologs. We observed that the PgCoV RBD bound human ACE2 at least as efficiently as the SARS-CoV-2 RBD, and that both RBDs bound pangolin ACE2 efficiently. We also observed a high level of variability in binding to closely related horseshoe-bat ACE2 orthologs consistent with the heterogeneity of their RBD-binding regions. However five consensus horseshoe-bat ACE2 residues enhanced ACE2 binding to the SARS-CoV-2 RBD and neutralization of SARS-CoV-2 pseudoviruses by an enzymatically inactive immunoadhesin form of human ACE2 (hACE2-NN-Fc). Two of these mutations impaired neutralization of SARS-CoV-1 pseudoviruses. An hACE2-NN-Fc variant bearing all five mutations neutralized both SARS-CoV-2 pseudovirus and infectious virus more efficiently than wild-type hACE2-NN-Fc. These data suggest that SARS-CoV-1 and -2 originate from distinct bat species, and identify a more potently neutralizing form of soluble ACE2.